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Abstract Colliding winds of massive binaries have long 
been considered as potential sites of non-thermal high- 
energy photon production. This is motivated by the de- 
tection of non-thermal spectra in the radio band, as well 
as by correlation studies of yet unidentified EGRET 7- 
ray sources with source populations appearing in star 
formation regions. 

This work re-considers the basic radiative processes 
and its properties that lead to high energy photon pro- 
duction in long-period massive star systems. We show 
that Klein-Nishina effects as well as the anisotropic na- 
ture of the inverse Compton scattering, the dominating 
leptonic emission process, likely yield spectral and vari- 
ability signatures in the 7-ray domain at or above the 
sensitivity of current or upcoming gamma ray instru- 
ments like GLAST-LAT. In addition to all relevant ra- 
diative losses, we include propagation (such as convec- 
tion in the stellar wind) as well as photon absorption 
effects, which a priori can not be neglected. 

The calculations are applied to WR 140 and WR 147, 
and predictions for their detectability in the 7-ray regime 
are provided. Physically similar specimen of their kind 
like WR 146, WR 137, WR 138, WR 112 and WR 125 
may be regarded as candidate sources at GeV energies 
for near-future 7-ray experiments. 

Finally, we discuss several aspects relevant for even- 
tually identifying this source class as a 7-ray emitting 
population. Thereby we utilize our findings on the ex- 
pected radiative behavior of typical colliding wind bina- 
ries in the 7-ray regime as well as its expected spatial 
distribution on the 7-ray sky. 
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1 Introduction 

By far the most convincing evidence for particle acceler- 
ation to relativistic energies mediated by the supersonic 
(terminal velocity Voo ^^1000-5000 km/s) winds of mas- 
sive (Af ~ 10-6 - -5MQ/yr), hot (T - 30000 - 50000K) 
stars comes from the observation of non-thermal radio 
emission (e.g. ]T|). This has been interpreted by syn- 
chrotron emission on the basis of the measured spectra 
(much steeper than the canonical value +0.6, F^, oc 
i/"' ) and high brightness temperatures of ~ 10^~^K, far 
exceeding ~ 10'*K expected from free-free emission from 
a steady-state isothermal radially symmetric wind [2]. 
Those particles have been suggested to be accelerated 
either in shocks caused by the instability of radiatively 
driven winds [5] , in the shocked wind collision region of 
multiple systems or in the termination shock [4^. Trig- 
gered by these observations it has been quickly realized 
that 7-ray emission should be expected as well, either 
through leptonic processes (inverse Compton scattering 
(IC) of the copious stellar UV photons (e.g. [S]), rel- 
ativistic bremsstrahlung (e.g. [5]) or hadronic interac- 
tions of co-accelerated ions with the dense wind material 
(e.g.[71[S]). This has established a plausible physical set- 
ting for massive star systems being putative 7-ray emit- 
ters. Indeed, positional coincidences of Wolf-Rayet (WR) 
stars with the population of so far unidentified EGRET 
sources have been found for 13 WR-binary systems [9l 
[TU1[TT| . Recently, the presence of non-thermal radio emis- 
sion has been linked to the binarity status of the stellar 
systems [12] , which supports the picture of particles be- 
ing predominantly accelerated at the forward and reverse 
shocks from the colliding supersonic winds from massive 
stars. 
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In this work we consider long-period binary systems 
as the most prospective 7-ray emitters detectable by the 
near-future experiment GLAST-LAT, and extend previ- 
ous theoretical work by including propagation, and an- 
isotropy and Klein-Nishina (KN) effects of the inverse 
Compton (IC) scattering process. 



2 The broadband SED of colliding v^rind systems 
of massive stars 

At present, the observationally established broadband 
spectral energy distribution (SED) of colliding wind bi- 
naries (CWBs) ranges from radio wavelengths to X-ray 
energies. At GHz frequencies, only one third of all ob- 
served WRs has been detected (e.g. [THIITI] ). Roughly 
40% of these show signatures of a contributing non-thermal 
component. These are predominantly binaries with peri- 
ods > 1 year, and this can be understood as an absorp- 
tion effect: The winds of massive stars are partially opti- 
cally thick at radio wavelengths. For typical 0-star wind 
parameters the = 1 surface lies roughly at ^ (1 — 2) x 
IO'^Rq at GHz frequencies. For WR stars these radii are 
even larger. Thus optically thin sight lines can only be 
found in systems with stellar separations of ^ 1O'^~'*R0. 
At IR wavelengths, massive binaries often show compo- 
nents in excess to the expected free-free emission (e.g. 
[15]). These persistent or variable/periodically/episodical- 
ly appearing IR excess fluxes have been interpreted as 
signatures of dust formation. The photospheric contin- 
uum emission dominates a broad range at optical to UV 
energies. Excess emission has also been observed at X- 
ray energies. Some systems show X-ray luminosities up 
to two orders of magnitude above the values suggested by 
the canonical relation Ljjf/Lboi ^ 10^'' obeyed by single 
early-type stars (e.g. [IS]). This excess emission has been 
attributed to the shock-heated (T ~ 10^~^K) plasma 
produced in the winds collision region. Indeed, recent 
Chandra observations of the close-by WR 147 system 
show extended emission peaking at a location which is 
in agreement with the expected wind collision region [17' . 
The observed phase-locked orbital variability in massive 
binaries in the soft X-ray as well as radio band is mainly 
a result of absorption in the dense stellar winds along 
the changing line of sight. 

In the following we shall concentrate on the non- 
thermal part of the continuum emission from CWBs. 



3 The model 

Various works have been devoted in the past to non- 
thermal emission components from CWBs (e.g. [5|[7lfT8l 
fT9ll20]. see also references in [21]) based on either lep- 
tonic (mostly IC) and/or hadronic (vr^-decay 7-rays) pro- 
cesses. They generally conclude on 7-ray luminosities in 
the range 10'^^ ' '^^erg/s to be expected from WR-binaries. 



Since the IC process has been shown to likely dominate 
the 7-ray production (e.g. [I9ll22]). this process is prone 
to determine the spectral appearance of WR-binaries at 
high energies. Here Klein-Nishina as well as anisotropy 
effects, neglected in past modeling of these systems, may 
provide valuable features that can be instrumental in 
identifying this source population at 7-ray energies. Fur- 
ther spectral imprints are expected from particle propa- 
gation within the extended colliding wind region. A sim- 
plified geometry of the system proofs sufhcient to high- 
light these effects. 

In the following we consider the sketch of a colliding 
wind region (CWR) that has been presented by e.g. Eich- 
ler & Usov (1993) T9] with the stagnation point defined 
by balancing the wind momenta with the assumption of 
spherical homogenous winds. Since in general Mwb. > 
Mob and Vqo^wb. ~ Voo,ob the shock distance to the OB- 
star, xoB, is much smaller than to the WR-star. We shall 
neglect here the interaction of the stellar radiation fields 
on the wind structure ( [23ll24j ): This effect is inherent to 
short-period binaries, and may influence the wind speed, 
thereby weakening the ram balance, shock strength and 
temperature. Thus our theoretical considerations here 
shall be restricted to long-period binaries. The stellar 
winds are permeated by magnetic fields originating from 
the surface of the massive star. Estimates for surface 
magnetic field strengths range from below Bg = 100 G 
(e.g. [25]) up to ~ lO-^G in WR-stars [2^. In the follow- 
ing we fix this value to a reasonable 100 G, unless stated 
otherwise, and use the well-developped magnetic rota- 
tor theory (e.g. [27]) to estimate the field strength Bq 
(in Gauss) at the CWR. Typically >mG or higher field 
strengths are expected, assumed to be constant through- 
out the emission region, at the CWR in long-period bi- 
naries. The shocked high-speed winds are creating a re- 
gion of hot gas that is separated by a contact discon- 
tinuity, and a forward and reverse shock follows. The 
gas flow velocity in this region away from the stagnation 
point will be some fraction of the wind velocity which 
we keep constant at V for simplicity. A simplification of 
the geometry from a bow-shaped to a cylinder-shaped 
collision region (with radius r perpendicular to the line- 
of-centres of the two stars, and given thickness) allows 
us to solve the relevant continuity equations analytically 
(see [28]). Here we consider diffusive shock acceleration 
(with acceleration rate a) out of a pool of thermal parti- 
cles, and take into account (continuous) radiative losses 
(synchrotron, IC, bremsstrahlung and Coulomb losses), 
(energy-independent) diffusion by introducing an escape 
time To and convection with speed V (set toV= 1/2wob 
if not noted otherwise) . The maximum particle energy is 
thereby determined self-consistently. 

At a distance > tq from the stagnation point con- 
vection along the post-shock flow will dominate over dif- 
fusion. At ro diffusion balances convection. Correspond- 
ingly, the emission region is devided into a region where 
acceleration/diffusion dominates, the " acceleration zone" , 
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Fig. 1 Energy loss rates in the acceleration zone due to 
IC scattering (thin solid lines) in the Thomson regime (TL) 
and KN regime (KN), synchrotron radiation (dashed-triple- 
dotted line) , relativistic electron-ion bremsstrahlung (dashed- 
dotted line) and Coulomb interactions (dashed line) in com- 
parison to the acceleration rate (thick solid line). Parameters 
are: bolometric OB-star luminosity Lboi.OB ~ 10"" Lq, target 
photon energy er ~ 10 eV, OB-star mass loss rate Mob = 
lO-''M0yr-\ WR-star mass loss rate Mwr = lOMoe, OB- 
star terminal velocity Uoo.ob ~ 4000km/s, binary separa- 
tion D — 10^'*cm, xoB ~ 0.24D, surface magnetic field 
Bs = lOOG, field strength at the CWR B « 0.5G, difltusion 
coefficient Ka = 1.6 ■ lO^^cm^s"^, escape time To ~ 1127 s, 
size of the acceleration region ro ~ 8.5 • 10^^ cm. See [5S] for 
details. 

and the outer region where convection dominates, the 
" convection zone" . The steady-state diffusion-loss equa- 
tion can be solved analytically provided suitable approx- 
imations for the KN cross section are applied (see [5H]). 
Fig. H] shows an example of all loss time scales for typical 
parameters of WR-binaries. The transition region, be- 
tween the Thomson and extreme KN range of the IC 
cross section, appears relevant for typical long-period 
WR-systems. A rigorous treatment of the Compton losses 
must therefore include KN effects. 

We find electron spectra with smooth roll-overs that 
cutoff at higher energies as compared to electron spec- 
tra that remain in the Thomson approximation through- 
out the whole particle energy range. In the convection 
region the particles loose energy radiatively as well as 
through expansion losses while in the post-shock flow. 
This leads to a dilution of the particle density as well 
as a deficit of high energy particles. The corresponding 
volume-integrated photon spectra softens (see Fig. 
depending on the relative size of the convection region 
with respect to the acceleration region. This itsself de- 
pends on the diffusion and convection time scales oper- 
ating in the CWR. 

With typically the photon to magnetic field energy 
density Uph/us ~ 52Lboi,OB,0/(-BG2^OB,0)^ > 1 (where 
iboi,OB,0 and xoB,Q are in units solar luminosity and ra- 



Fig. 2 IC spectra from the acceleration region for D = 
10^'' cm and &l ~ 0° (lower solid line), 30° (lower dashed line), 
60° (dotted line), 90° (dashed-dotted line), 120° (dashed- 
triple dotted line) 150° (long dashed line) and 180° (upper 
solid line). For 9^ ~ 30° the total volume- integrated (i.e. ac- 
celeration plus convection zone) IC spectrum is also shown 
(upper dashed line). All other parameters are the same as in 

Fig.m 

dius, respectively) IC scattering turns out to be the most 
important radiative loss channel for relativistic electrons 
in CWBs in most cases [T9lf22] . Because the stellar target 
photons for IC scattering arrive at the collision region 
from a prefered direction, the full angular dependence 
of the cross section has to be taken into account (e.g. 
|29j). This leads to anisotropy effects like the emitted 
flux and cutoff energy being dependent on the sight line 
into the wind (see Fig. [2]). The calculations are carried 
out using the monochromatic approximation for the pho- 
tospheric target field and an isotropic particle distribu- 
tion in the emission region. The latter is maintained by 
pitch angle scattering at magnetic field inhomogeneities. 
Fig. [2] shows the maximum flux level reached at phases 
where the WR-star is behind the OB-star along the sight 
line. Furthermore, large scattering angles tend to pro- 
duce more energetic photons than low angles. This can 
affect the cutoff photon energy by several orders of mag- 
nitude for a given value of the electron energy cutoff (see 
Sect. 4.2). 

Details are described in [5B]. 



4 Application to archetypal systems 

4.1 WR 140 

The archetypical WR-binary system WR 140 (WC 7pd 
-I- 04-5 V), located in the Cygnus constellation at a dis- 
tance of ~ 1.85 kpc, is one of the most detailled studied 
specimen of its kind. Its long ('^ 8 years) period and ex- 
treme eccentricity (e « 0.88) makes it a diverse system 
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to study its non-thermal behavior in the radio band. For 
the first time, Dougherty et al. (2005) |30| resolved the 
bow-shaped arcs of the emission region at 8 epochs of 8.4 
GHz VLBA observations. We used the synchrotron spec- 
tra and system parameters as published in [30j to predict 
its non-thermal high energy emission at various phases. 
With these values the CWR is located at a distance xqb 
of 0.32 times the stellar separation. For a surface mag- 
netic field of 100 G roughly equipartition values for its 
field strength at the CWR location follow. With a rela- 
tivistic particle injection energy of ^ 10~^'"~^% of the 
kinetic wind energy, the synchrotron fluxes at the con- 
sidered phases (0.2, 0.67, 0.8, 0.95) could be reproduced. 
From radio observational grounds, relativistic electrons 
at least up to ^10-100 MeV do exist. This limits the dif- 
fusion coefficient n to sufficiently low values which will 
allow the acceleration rate to overcome the Coulomb loss 
rate at low energies. A further constraint of the diffusion 
coefficient is provided by the spectral shape of the non- 
thermal radio component, if the shock compression ratio 
and convection velocity are known. On the other side, k 
is limited by Bohm diffusion. For k — 2 ■ lO^^cm^s"^ up 
to ~ lO^MeV-electrons are expected at least at apastron, 
where radiative losses are sufficiently small not to affect 
the cutoff particle energy (see Fig. [3]) . The self-consistent 
determination of the maximum particle energy allows 
sound predictions at the highest energies. Note that KN 
effects will alter the emitting electron spectrum above 
~ lO^'MeV, as indicated in Fig. El This is in contrast to 
phases close to periastron. The intense stellar radiation 
field there limits the electron spectrum to ~ lOOMeV 
(see Fig.[3|). As a consequence, the corresponding photon 
spectrum from IC scattering cuts off already in the soft 
7-ray band (see Fig.[4|), and hadronically produced pho- 
tons may dominate at (sub-)GeV energies, albeit with a 
possibly undetectable flux level even for the more sensi- 
tive near- future experiments. EGRET observations, that 
were carried out rather around periastron, therefore did 
not lead to detections. Mainly as a result of including 
also hadronic ion-ion interactions into the calculations 
for TT^ production, the corresponding Tr'^-decay 7-ray flux 
in Pittard & Dougherty [31] extends to higher energies 
than considered here. The association of WR 140 to the 
unidentified EGRET source 3EG J2022-f 4317 appears to 
be vague: With WR 140 being located ~ 0.67° away from 
the nominal position of 3EG J2022-I-4317, and barely 
consistent with the 99% source location uncertainty con- 
tour a conclusive identification seems farfetched. It is 
due to the extreme eccentricity of this system, leading 
to significant changes in the stellar radiation field den- 
sity at the shock location, that causes a blurring of the 
phase-locked flux variations expected otherwise from the 
anisotropy of the dominating IC scattering process. In- 
deed, orbital flux variations in this system are reduced 
to a factor ~ 2 — 3 (see Fig. |4|). In general, relativistic 
bremsstrahlung radiation lies always below the IC emis- 
sion level in these systems. When compared to the ex- 
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Fig. 3 Steady-state electron spectra for WR 140 at orbital 
phases 0.955, 0.2, 0.671 and 0.8. Flux is in arbitrary units. 
See text for parameters. 

pected IC flux, Tr'^-decay 7-ray production is small, even 
if the total wind energy is transformed into relativis- 
tic protons/ions. This makes WR-binary systems rather 
unpromising putative neutrino and cosmic ray sources 
within our model. 

Fig. [4] implies that WR 140 may be detectable with 
GeV-instruments like GLAST-LAT even at individually 
selected phases if the electrons reach sufficient high ener- 
gies, while INTEGRAL requires Msec exposures for any 
detection above 1 MeV. The expected extent of the pho- 
ton spectrum to ^ 10 — 100 GeV at phases where the 
binary separation is large, may even allow low-energy 
threshold Imaging Atmospheric Cherenkov Telescopes 
(lACTs) to gain important information on the cutoff 
of WR 140's spectrum. This may allow to further con- 
strain the diffusion coefficient. Note that photon absorp- 
tion due to 77 pair production can not a priori be ne- 
glected in colliding wind systems. E.g. for WR 140 the 
optical depth T~^^{r = ro,E « lOOGeV)^ 1 at phase 
0.67, T^^{r ro, S « lOOGeV)- 3 at phase 0.95. 

4.2 WR 147 

Due to its proximity WR 147 (WN8h -I- B0.5V) is one 
of the few binary systems where high resolution radio 
observations lead to resolving this system into its com- 
ponents (Williams et al. 1997). The observed radio mor- 
phology supports a binary separation of ~ 417 AU for a 
source distance of 650 pc. We used the system parame- 
ters as published in Setia Gunawan et al. (2001) to- 
gether with the observed synchrotron spectrum to model 
its expected spectral behaviour at high energies. Since 
WR 147's inclination i nor eccentricity e are known, we 
assumed for our modeling e = and i — 90°. The ob- 
served synchrotron spectrum could be reproduced within 
its observational uncertainties if ^ 0.15% of the OB- wind 
kinetic energy is transformed into relativistic electrons. 
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Fig. 4 IC spectra for WR 140 at phases 0.955, 0.2, 0.671 
and 0.8 from electron spectra as shown in FigO The spec- 
tral changes from 7-ray a bso rption are not shown here. The 
EGRET 2(7 upper limit [22] is based on observations that 
correspond to a superposition of orbital states rather deter- 
mined by periastron phase, and therefore applies rather to 
phase 0.955. 



and assuming a surface magnetic field of 30 G (translat- 
ing into 25 mG fields at the CWR). Due to its huge bi- 
nary separation the emitting particle spectra are not lim- 
ited by IC losses but rather by the size of the acceleration 
region. Fig. [5] presents the expected orbital IC flux vari- 
ations due to the anisotropic nature of the IC process for 
the above described parameter set. Here the maximum 
flux and photon cutoff energy, expected when the WR- 
star is behind the OB-star along the sight line, are more 
than one order of magntitude higher than their minimum 
values. Absorption of > 50 GeV photons turns out neg- 
ligible for this system, provided its eccentricity is small. 
The chance to detect WR 147 with future high sensitive 
instruments is somewhat better than for WR 140. While 
LAT will probe the MeV-GeV range of its spectrum, 
only the low threshold (< 100 GeV) experiments in the 
northern hemisphere among the ground-based lACTs 
may have a distinct chance to detect WR 147 at its high- 
est energies. 



4.3 Further candidates ? 

Our work provides ample evidence for the archetypal 
long-period WR-binaries WR 140 and WR 147 being 
plausible candidates for detection with near future 7-ray 
instruments. Both provide target photon fields at the 
CWR which is sufficiently dense to produce a IG scat- 
tered photon flux above the respective instrument sen- 
sitivity, but at the same time does not act as a thick 
absorber for those photons, nor prevents electron accel- 
eration to high enough energies by radiative losses. A 
sensible search for further 7-ray source candidates may 
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Fig. 5 IC spectra of WR 147 for orbital phases 0, 0.25, 0.5 
and 0.75. See text for parameters. 77 pair production absorbs 
not more than <0.3% (> 50 GeV) and <18% (> 100 GeV) of 
the produced flux at orbital phases 0.25 and 0.5, respectively 
(not shown in figure). No absorption takes place at phase 0. 



therefore start with scanning the WR-binary population 
for physically similar systems. A (not complete) list en- 
compasses WR 146, WR 137, WR 112 and WR 125, 
with WR 112 being the most distant (~ 4 kpc) of them. 
From all of them non-thermal radio components have 
been observed at least occasionally. One (WR 137) has 
been found positionally coincident with an unidentified 
EGRET source (3EG J2016-h3657: 0). 



5 Perspectives for the gamma-ray band 

The identification of a specific colliding wind system or 
the population thereof as 7-ray emitters requires first: a 
convincing positional association, and second: its physi- 
cal relation with the detected 7-ray source. 

The most up-to-date catalog of WR-stars in the Milky 
Way lists more than 220 sources with a detected bi- 
nary frequency of ~ 40 — 50% in the solar neighbor- 
hood [33]. Its spatial distribution reflects the spiral arm 
structure of our Galaxy with a slight asymmetry with re- 
spect to the galactic plane (^33^ and references therein). 
64 WR- (systems) are found in the solar neighbourhood 
(< 2 kpc), more than 25 have been detected in the im- 
mediate vicinity (< 30pc) of the galactic center. The 
total number of WR- (systems) in the Galaxy has been 
estimated to up to ~ 8000 [34], total OB-star number 
counts in our Galaxy may reach values as high as 60000 
[35j . This is orders of magnitude larger than the number 
of all, presumably galactic, but still unidentified 7-ray 
sources detected to date. 

Intriguing positional coincidences of massive star pop- 
ulations, among them colliding wind systems, with uniden- 
tified EGRET sources have been found in the past (e.g. 
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[9|ll0j). however no individual binary system could un- 
ambigeously be identified as a high energy 7-ray emitter 
so far. This is most likely a statistical issue. EGRET's 
large location error box, typically ^ 0.25° for a strong 
point source at high latitudes, and even larger for sources 
in the galactic plane, allows to find many plausible coun- 
terparts within the positional uncertainty of the 7-ray 
source. Furthermore, the strong diffuse 7-ray background 
at low galactic latitudes makes the detection of any non- 
periodic galactic source close to the instrument sensitiv- 
ity even more challenging. The capabilities of EGRET's 
successor, GLAST-LATQ, will improve here. Its expected 
ability of both detecting weak signals (> 3 • 10~^cm~^ 
s""'^) and at the same time locating an object in the 7-ray 
sky within ~ 0.5° — 5', solves the confusion problem at 
least among the galactic bright EGRET sources. 

How will a typical CWS appear at 7-ray energies? 
What are its observable characteristics that may discrim- 
inate it from the wealth of other plausible 7-ray source 
classes? 

To begin with, the available energy kinetic wind en- 
ergy = Mv'^/2 provides a firm upper limit on the 
photon production rate from each process, and thus on 
the expected flux level. Inserting typical mass loss and 
wind parameters one finds the kinetic energy of massive 
star winds to be of order 1% of the bolometric radiative 
energy output, or typically ~ lO'^^erg/s. Taking into ac- 
count the efficiency of the particle acceleration process, 
much smaller values will rather be the rule. 

Furthermore, colliding wind regions will appear not 
extended, but point-like on the 7-ray sky for both present 
and near future satellite and ground based instruments. 

Our study implies that detectable phase-locked or- 
bital variations in the 7-ray band occur, however, they 
may possibly be blurred by the geometry of the system 
(in particular due to eccentricity and inclination). Flux 
variations at high energies may additionally be expected 
from modulations of the geometry and size of the emit- 
ting region, modulations of the stellar target photon field 
density, inhomogeneities and/or shocks in the wind out- 
flow, etc. The time scale of these variations is determined 
by the respective physical cause. 

Spectral hints on the expected low- energy IC compo- 
nent of CWBs may be deduced from their visible syn- 
chrotron spectrum. Spectral indices at GHz frequencies 
in the range ar=0.3...-1.3 has been detected so far from 
the current limited sample of non-thermal emitters among 
the WR-binary systems. This may hint towards photon 
spectra with photon index aph < 2.5 far below the cutoff 
energy. At GeV-energies KN, anisotropy and propaga- 
tion effects may have an impact on the spectral shape. 

With these characteristics at hand, a population study 
may be instrumental to finally unveil the class of massive 
star binary systems as high energy emitters. 



^ http:/ /www-glast. slac.stanford.edu/software/IS/glast_lat_i 
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